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We show that time-reversal invariant superconductors in d = 2 (d — 3) dimensions can support 
topologically stable Fermi points (lines), characterized by an integer topological charge. Combining 
this with the momentum space symmetries present, we prove analogs of the fermion doubling theo- 
rem: for d — 2 lattice models admitting a spin x electron-hole structure, the number of Fermi points 
is a multiple of four, while for d — 3, Fermi lines come in pairs. We show two implications of our 
findings for topological superconductors in d — 3: first, we relate the bulk topological invariant to a 
topological number for the surface Fermi points in the form of an index theorem. Second, we show 
that the existence of topologically stable Fermi lines results in extended gapless regions in a generic 
topological superconductor phase diagram. 
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I. INTRODUCTION 

Time-reversal (T) invariant superconductors have at- 
tracted a lot of attention lately as they provide possi- 
bility to realize topologically nontrivial gapped phases 
of matter»ir— While topologically nontrivial T-invariant 
gapped phases of 3 He were discussed already back in 
the eighties^— (these phases also appear in the re- 
cent constructions^ see also Ref. Qjl), the key obser- 
vation highlighting the general role of T and electron- 
hole (EH) symmetries in the topological classification of 
gapped T-invariant superconductors has been made only 
recently: Schnyder et a& has noted that combining T 
and EH symmetries, one finds that the Bogoliubov dc 
Genncs Hamiltonian 

(H(k) A(k) \ m 

obeys a chiral symmetry: CfiBdcC 1 = —HbcIGi with 
C = C f unitary. In Eq. <JTJ) , the N x N matrices HQs.) and 
A(k) are the single particle Hamiltonian and the pairing 
potential. There exists a basis 1 - in which the chirality 
operator is C = diag(ljV) — 1jv)i an d the Hamiltonian 
takes the block offdiagonal form, 

?4d G ( k ) = (^tjk) ^o k) ) ■ ( 2 ) 

A gapped system, as the determinant of D(k) is nonva- 
nishing for all k, realizes a mapping from the momen- 
tum space to the space GL(N,C) of invertible N x N 
matrices. In d — 3, such mappings fall into distinct ho- 
motopy classes, which directly translates to the possibil- 
ity for topologically nontrivial three dimensional gapped 
systems, i.e., systems which are not dcformable into a k- 
independent gapped model without a gap closing some- 
where throughout the deformation^ Such topologically 
nontrivial gapped superconductors are called topological 
superconductors. (Topological superconductors also ex- 
ist in d = 1, 2 fir— ^ but the ones considered in this work 
are three dimensional.) The homotopy class of the sys- 



tem is characterized by an integer bulk topological invari- 
ant constructed such that its nonzero value signals a 
topological superconductor. 

Besides their value in analyzing gapped phases, ho- 
motopy arguments prove to be useful in studying sys- 
tems without a gap»i^— In such systems, homotopy ar- 
guments can be used to determine the possible dimen- 
sionalities of topologically stable Fermi surfaces, defined 
as manifolds of zero energy that cannot be destroyed by a 
small deformation of the Hamiltonian. A Fermi surface is 
stable if a nontrivial topological charge can be assigned 
to it; the possible values of the topological charge also 
result from the homotopy considerations. 

As we have illustrated, the block-offdiagonal struc- 
ture in Eq. @ provides key insights into the topology 
of gapped phases of T-invariant superconductors. In 
this work we would like to study what insights it could 
provide regarding the nature of gapless phases of these 
systems. One motivation for this, naturally, lies in the 
characterization of the possible nodal structures of sys- 
tems such as high-T c materials, or noncentrosymmetric 
superconductors^ The reason, however, why this ques- 
tion is especially timely is given by the relevance of gap- 
less systems in the context of topological superconduc- 
tors: first, as v does not change under non-gap closing de- 
formations of the theory, the system must become gapless 
across a topological phase transitio n 14 ' 16 (a transition be- 
tween gapped phases with different v). Second, the sur- 
face of a topological superconductor can also be thought 
of as a two dimensional gapless system: it is known to 
support topologically protected gapless modes^i^ robust 
against arbitrary deformations (provided that these do 
not close the bulk gap). 

With these motivations in mind, our first objective is 
to present a homotopy argument leading to the topolog- 
ical characterization of Fermi surfaces of d = 2, 3 dimen- 
sional T-invariant superconductors. This is carried out 
in Sec. [Til In Sec. IIIII we deduce the consequences of our 
findings for two dimensional systems, including the sur- 
face of a topological superconductor. The consequences 
in three dimensions are discussed in Sec. IIV1 together 
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with the implications regarding topological phase transi- 
tions. Our conclusions are summarized in Sec. [V] 

II. FERMI SURFACES AND TOPOLOGICAL 
CHARGE 

Due to det H^ dG (k) = — | dct£>(k)| 2 , the vanishing of 
the gap at a point ko on a Fermi surface is in one-to-one 
correspondence with the vanishing of the determinant of 
D(k) at this point. We first assume that ko is not a 
time- reversal invariant momentum (TRIM), i.e., there is 
no such reciprocal lattice vector G that ko = — ko + G. 
Then, in the vicinity of ko, the matrix -D(k) need not sat- 
isfy any restrictions; all restrictions stemming from EH 
and T symmetries relate^ -D(k) to D(— k). Therefore, in 
the vicinity of k the only restriction on T-L' BdG (k) is chi- 
ral symmetry. The properties of Fermi surfaces in d = 2 
systems with this restriction were studied by Wen and 
ZeeJ£ Their considerations can be formulated in a com- 
pact way which also includes the d = 3 case as follows 
(see also Refe. [IMfi!): for a Fermi surface of dimension 
d — p, in the p dimensions transverse to it, det -D(k) ^ 
in the vicinity of the surface. Then, embracing the point 
ko on the Fermi surface with a small p — 1 sphere S p ~ l 
in the transverse dimensions, we have a mapping from 
S' p_1 to the space GL(N, C). If this represents a nontriv- 
ial class in the homotopy group tt p _i(GL(N, C)), there is 
no small deformation of -D(k) which would make the map 
contractible to the inside of S' p_1 , thus would remove the 
point of det D(k) — within S^ -1 , i.e., the Fermi surface 
is stable. For p = 1, 2, 3 only for p = 2 is there a non- 
trivial homotopy group, tt 1 (GL(N,C)) = Z. This leads 
to our central result: T-invariant superconductors can 
support d — 2 dimensional topologically stable Fermi sur- 
faces, which can be assigned a topological charge taking 
values in the additive group of integers. 

Fermi surfaces in d = 3 (d = 2) are thus Fermi lines 
(Fermi points). Their dimensionality can be also under- 
stood by noting that a Fermi surface is a manifold of 
dctD(k) = 0, i.e., it is defined by two real equations for 
k.— The topological charge can be found by noting that 
the sphere in the transverse space is S , i.e., it is a cir- 
cle. The topological invariant q characterizing mappings 
from the circle to GL(N, C) is simply given by the wind- 
ing number of the phase of det£>(k). The topological 
invariance of this charge means that q docs not change 
under deformations of the circle into any loop around 
the Fermi surface, nor under deformations of D(k), pro- 
vided that the loop does not cross a Fermi surface during 
these changes. Following Ref. [H (see also Refs. Il6ll20| ). 
we can express this integer charge directly through the 
Hamiltonian as 

q= 4^f Tl Wado)' 1 W BdG ] dl, (3) 

where the integral is for a loop around the Fermi surface. 

Due to its dimensionality, and due to the fact that q is 
the winding of the phase, a Fermi surface can be viewed 



as a momentum space vortex in detD(k) as an "order 
parameter" . Henceforth, we will often use this vortex 
picture; we emphasize that in this work, vortices, antivor- 
tices always refer to these momentum space objects and 
not to real space defects in the superconducting order pa- 
rameter. Under a deformation of D(k), locate creation 
of Fermi surfaces proceeds through formation of vortex 
loops in three dimensions, and vortex-antivortcx pairs 
in two dimensions. Processes in the opposite direction 
are also possible: a large enough deformation can shrink 
and remove vortex loops, or annihilate vortex-antivortex 
pairs; this explains why the notion of topological stabil- 
ity provides robustness only against small perturbations. 
For lattice systems, instead of referring to the way the 
Fermi surfaces are created, a more general statement can 
be formulated, similarly to Ref. I22I: assume that D(k) is 
smooth and periodic in momentum space. In d = 2, the 
net winding is then zero in the reciprocal unit cell and is 
equal to the sum of the topological charges of the Fermi 
points. In d = 3, a similar statement can be made for 
two dimensional sections of the reciprocal unit cell paral- 
lel to its faces, with the intersections of Fermi lines with 
a section playing the role of Fermi points. 

At this stage, it is worthwhile to compare our result 
to earlier findings for single band systems * 16 ' 7 ' 23 While 
these works also find d — 2 dimensional Fermi surfaces, 
Ref. [I?] finds a topological charge with Z2 values in- 
stead of Z values. (Ref. [l6| found a Z charge for a spin- 
degenerate, single band system.) For comparison, we 
consider the two-dimensional example of high-T c mate- 
rials from Ref. The single particle Hamiltonian and 
the pairing potential are 

H(k) = a(k 2 - k 2 F ) + g(k)<7, (4) 
A(k) = iA (fc^-fc2) ( T 2 +id(k)cr, (5) 

where Oj are the Pauli matrices in spin-space, and the 
band parameter a, the singlet pairing strength Ao, and 
the Fermi momentum hp are positive constants. The vec- 
tors g(k) and d(k), corresponding to spin-orbit coupling 
and the triplet pairing strength, respectively, are real odd 
functions of k. For g(k) = d(k) = 0, there arc four Fermi 
points at ko = (±fcp, ±kp)/y/2, which split into pairs of 
Fermi points once g(k), d(k) become nonzero. Ref. [i~7l 
finds that the Fermi points for g(k) = d(k) = have 
zero charge and they split up into pairs with unit charge, 
in agreement with the Z2 sum rule. As a Fermi point 
with zero charge is not topologically stable, an alterna- 
tive scenario seems also possible, in which the system 
becomes gapped due to a small deformation. Apply- 
ing our considerations shows that this is not the case. 
After a suitable transformation diagonalizing the chi- 
rality operator, the offdiagonal block of T-L' BdG ()i) be- 
comes D(k) = A(k) + H(k)a 2 . For g(k) = d(k) = 0, 
det .D(k) = £(k) 2 where the winding of the phase of £(k) 
around each Fermi point is of unit magnitude, hence for 
each Fermi point |g| = 2. Under a small deformation, the 
splitting into two Fermi points with \q\ = 1 is thus the 
generic behavior. The example here suggests Ref. [l?] to 
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be consistent with our findings upon identifying the Z2 
charge as q mod 2. 

III. CONSEQUENCES IN TWO DIMENSIONAL 
MODELS 

A. Fermion quadrupling in lattice models 

After discussing aspects for which the restrictions re- 
lating -D(k) to D(— k) arc unimportant, we turn to study 
the implications of these restrictions. Time-reversal in- 
variant superconductors can be classified depending on 
symmetries under spin-rotation: they are said to belong 
to class Dili if there is no spin rotation symmetry, to class 
AIII if 03 is conserved and to class CI if there is full spin- 
rotation symmetry^ - — The restrictions on -D(k) are as 
follows^: 

D(k) = -D T (-k) (class Dili), (6a) 
A^k) = Dj(-k) (class AIII), (6b) 
Ax(k) = D s (k) = D^(-k) (class CI). (6c) 

For class AIII, a chiral symmetry C = diag(ljy/27 — ljv/2) 
applies for each spin separately. From these relations 
we see that if ko is on a Fermi surface, so is — ko. Re- 
garding the further consequences of Eqs. ©, we first fo- 
cus on the two dimensional case, i.e., the Fermi surfaces 
are Fermi points. We consider the generic situation that 
Fermi points have \q\ = 1 (not counting spin in class CI); 
a Fermi point with a higher \q\ splits into Fermi points 
with \q\ = 1 upon a slight deformation of D(k). We start 
with class Dili and CI. Let us first assume that ko is 
not a TRIM. Taking a loop around ko, and the momen- 
tum reversed of this loop around — ko, calculating the 
phase winding of dctD(k) on these loops, we immedi- 
ately find that q of a Fermi point at ko and at — ko is 
the same. What happens if ko is a TRIM? For class CI, 
for each point on a circle centered at a TRIM, the de- 
terminant of -D(k) at the antipodal points of the circle 
agree, leading to an even winding number. Thus, if there 
is a Fermi point at a TRIM, a deformation of D(k) splits 
it into Fermi points off the TRIM, leading us back to a 
case satisfying our original assumption. The same rea- 
soning can be applied for class Dili, provided that D(k) 
is even dimensional. This is guaranteed if % BdG admits 
a spin x EH structure, since it implies that its dimension 
is a multiple of four. For class AIII, if there is a Fermi 
point at ko for spin up, there is an other at — ko for spin 
down. 

If the conditions for zero net winding hold, each vor- 
tex should have its antivortex (for class AIII, for each 
spin separately). It follows from our preceding discus- 
sion that the momentum reversed pairs cannot be also 
the antivortex pairs. Starting from a Fermi point at ko, 
we thus immediately obtain three other. That is, the 
combination of Eqs. ([5]) with our findings regarding the 
topology of Fermi points has led us to a result analogous 



to the fermion doubling theorem^ the number of Fermi 
points is a multiple of four; for class CI each Fermi point 
is spin degenerate (i.e., we have a multiple of eight in- 
cluding spin), for class AIII, the number of Fermi points 
is a multiple of two for each spin. We note that Ref. 
conjectures a similar result, stating that the number of 
Fermi points is an even multiple of the minimal number 
of Dirac cones in a given symmetry class. This agrees 
with our findings for classes AIII and CI, but it gives 
the number of Fermi points to be a multiple of two for 
class DHL With our proof at hand we see that this is not 
possible in models with a spinxEH structure. 

B. Surface states of topological superconductors 

A two dimensional surface separating topological su- 
perconductors with different bulk topological invariants 
v+ and is known to support topologically protected 
Fermi points*^ (The theory for the low energy exci- 
tations on the surface defines a two dimensional T- 
invariant superconductor model, hence its gapless modes 
form Fermi points.) Considerations based on three di- 
mensional Dirac models^ suggest the number of surface 
Fermi points to be | v + — \ . This equality cannot hold 
in all cases: while v+ — v- is a topological invariant, 
i.e., it cannot change under deformations (which do not 
close the gap), the number of surface Fermi points is 
not, as such deformations are allowed to add or remove 
vortex-antivortex pairs in the surface theory. Note, how- 
ever, that similarly to the violation of the fermion dou- 
bling theorem on the surface of a topological insulator^ 
the fermion quadrupling theorem can be violated in the 
above domain wall by allowing the Fermi points to have 
a nonzero net winding number <7 net (i.e., it is possible 
to have vortices without antivortex pairs). This num- 
ber, as it is quantized in integer values, also cannot 
change by continuous deformations of the theory: it is 
also a topological invariant. (Stated in another way, q ne t 
measures the number of unpaired vortices; this clearly 
cannot change by adding or removing vortex-antivortex 
pairs.) It is tempting to look for an "index theorem" 
relating these topological invariants. For this relation, 
we can use the results from the Dirac models as a ref- 
erence. For a domain wall in the x — y plane with 
v{z — > ±00) = v±, the \v+ — Fermi points of the 
surface theory are characterized^ by a D(k) equivalent 
to (fci + Tjik2)ti v v _\, where 77 = sgn(z/ + — vJ). Im- 
portantly, all Fermi points are unpaired: they all have 
identical unitary winding q = sgn(z/ + — The net 

winding is thus 

feet =v+-v-, (7) 

providing us with the index theorem we were looking for. 
This equation, due to the topological invariance of both 
sides, holds irrespective to the possible changes in the to- 
tal number of surface Fermi points. For this total num- 
ber, assuming again that all Fermi points have \q\ = 1, 
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| v + — V- 1 sets the minimum: once all the Fermi points 
are unpaired, their number cannot be further reduced 
by pair annihilation. This also immediately explains the 
topological protection of the gapless surface states. 

As a special case, Eq. (JT]) includes the surface of a 
topological superconductor, i.e, an interface between a 
topological superconductor (y = v + ) and the vacuum 
(z/_ = 0). In this case, the significance of our index theo- 
rem is that it gives a direct measure of the bulk topolog- 
ical invariant of the superconductor in terms of a topo- 
logical number for the surface degrees of freedom. 

It is worthwhile to go through the possible values of 
(/net for the various symmetry classes in Eq. ([6]). Our 
reasoning in Sec. IIII Al shows that g ne t is always even in 
class CI (not counting spin). This is in accord with the 
observation 6 - that v can take only even values in this sym- 
metry class. For classes Dili and AIII (considering only 
one spin in the latter case), v can take any integer value. 
This is again in accord with our previous considerations 
if we note that for class Dili, the surface theory model 
is not constrained to have a spinxEH structure. 

Finally, our index theorem can be also viewed as com- 
pleting an analogy between (T-breaking) chiral super- 
conductors in two dimensions, and T-invariant topologi- 
cal superconductors in three dimensions. Regarding bulk 
properties, the two classes of systems are analogous in 
that both are characterized by a bulk topological invari- 
ant v taking values in Z. For a domain wall in chiral 
superconductors, there is an index theorem^ equating 
v+ — V- with the number of unpaired edge states (i.e., 
the difference between the number of left moving and 
right moving modes in the domain wall). Our index the- 
orem with the number of unpaired Fermi points on its 
left hand side thus nicely completes the analogy for the 
boundary degrees of freedom. 



IV. CONSEQUENCES IN THREE DIMENSIONS 

A. Fermi lines come in pairs 

In three dimensions, a d — 2 dimensional Fermi surface 
is a Fermi line. We now prove our next result: the conse- 
quence of Eqs. © is that Fermi lines come in pairs. While 
for class AIII, this is simply due to a Fermi line with spin 
up having a pair with spin down, classes CI and Dili de- 
serve some discussion. We consider a lattice system with 
reciprocal unit cell defined by a.ih i: a; G (—1/2, 1/2). 
(Here bi are primitive reciprocal lattice vectors). The 
half cell bounded by the plane V at — and the face 
V' with »3 = 1/2 is mapped to the other half of the cell 
upon k — > k. Due to Eqs. (O, a Fermi line located 
entirely in one of the half cells (not intersecting V nor 
V') has a pair in the other half cell. The situation is 
not so clear for Fermi lines intersecting V and/or V'\ to 
prove our assertion we now show that in the generic situ- 
ation a Fermi line cannot be its own pair under k — > — k. 
First note the following facts: the TRIM of our recipro- 



cal unit cell are on V ', V' . The behavior of V , V' under 
k — > k is that of the two dimensional case, thus the 
intersection points with Fermi lines have the earlier dis- 
cussed properties of Fermi points in two dimensions. In 
particular, generically there are no intersection points at 
TRIM. Now assume that there is a Fermi line which is its 
own pair, which intersects V and V in the set of points 
{k,-} U {k'\ U {-k,-} U { K\ (k, G V. kj G V). Con- 
sider the pair of points ki, — ki with winding number 
q. Take two parametrizations of the line, a(i) and b(i), 
a(0) = ki, b(0) = -k x , a(t) = -b(t). Since the Fermi 
line is a vortex, for a given parametrization, Fermi points 
on V , V with like charges are pierced from the same di- 
rection (always from above or always from below). As ki 
and ki have the same charges and since b pierces — ki 
from the opposite direction to how a pierces ki , a and b 
parametrize the line with opposite orientation. Thus the 
two parametrizations must run into each other at some 
to, where a(io) = h(to) + G, which implies that a(io) is 
a TRIM, which does not happen in the generic situation. 
This thus proves that lines come in pairs also for classes 
Dili and CI. While for class CI, it would also possible, as 
for class AIII, to define pairs in terms of spin, as we saw, 
more can be said: there arc spin degenerate pairs of lines. 
Note that in our proof we have not used the criterion of 
zero net winding. This however does not lead to further 
doubling of the number of lines, as it can be satisfied for 
example with pair of Fermi loops. 



B. Phase diagram of topological superconductors 

A topological phase transition, i.e., a transition be- 
tween gapped phases with different bulk topological in- 
variant v proceeds through a gapless phase. 14 ' 16 If the 
gapless phase can support topologically stable Fermi sur- 
faces, it is expected to be a stable region in parame- 
ter spaced small variations of the system parameters 
should not lead out of it. In earlier models of topolog- 
ical superconductors j 6 - - — however, topologically distinct 
gapped phases are directly adjacent to each other in pa- 
rameter space; a slight deformation of a gapless system 
on the phase boundary between them immediately leads 
to one of the gapped phases. The gapless phase is thus 
unstable. As a final application, we show how the exis- 
tence of topologically stable Fermi lines can modify this 
situation. We consider the class CI model of Schnyder et 
al,— constructed for the diamond lattice. It is defined by 
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Figure 1: The phase diagram of the model in Eq. (|10|) . The 
dashed lines are the phase boundaries for /i = [i.e., for the 
model in Eq. ([8)l] running at fi s = ±6£'. The solid lines are 
the numerically obtained phase boundaries for t = 4, A = 2, 
= 0.5; they enclose a gapless region in parameter space. 
With dots, we also displayed /j, s — ±6t' ± fi, defining the 
phase boundaries according to Eq. (jT2j ■ 



0k = 4i' COS — f COS — — COS — I + u s , 
2 V 2 2 J P ' 

where the real numbers i, t' parametrize the nearest 
neighbor and next-nearest neighbor hopping amplitudes, 
A is a real pairing strength, and ll s is a staggered poten- 
tial. The phase diagram is shown in Fig. [TJ The system 
is gapless on the lines ll s = ±6i', with the gap closing oc- 
curring at momenta K\ t ± = 2tt ( ±1/3, 1, ) , K% t ± = 
2ir ( 1, ±1/3, ). Gapless systems belong to an unsta- 
ble region in parameter space: a slight change for exam- 
ple in ll s immediately opens a gap. This can be traced 
back to the fact that the model satisfies 

D(k) =T 2 D*(k)T 2 , (9) 

where Tj are Pauli matrices in the grading of Eq. ©. 
Consequently, detD(k) = A 2 . + |$ k | 2 + 6|, thus the 
space of invertible D(k) matrices is R 4 \ {0}. As 
7T p _i(]R 4 \ {0}) = (p = 1,2,3), the model is indeed not 
expected to have stable gapless phases. With Eq. © 
in addition to Eq. (|6cl) . the model has more symmetries 
than required for belonging to class CI (i.e., it is "not 
generic"). The symmetry §§§ can be removed by intro- 
ducing a nonzero chemical potential ll. (The ll = case 
corresponds to half-filling.) This leads to 

D'(k) = D{k) + i^ 2 , (10) 

which, due to dot D(k) = (A k + i/i) 2 + |$ k | 2 + O 2 ., mod- 
ifies the space of invertible matrices. To see how the 
nature of the gapless phase is modified, we assume that 
if, li Si //is small compared to t, A, and expand -D'(k) 
around the points Kj t ±, resulting in 

D' h ±(p) = (pi + tM)l2 + iriP2 + ir 2 Pz - «T 3 m ± , (11) 



where p is a scaled momentum parametrizing the devia- 
tion from Kj^±, and rh± = fx s ± 6t'. The condition for a 
gapless phase is now 

H 2 -m% =p 2 2 +pl pi=0, (12) 

describing the appearance of Fermi lines through the 
transition. In this more generic, ll ^ case, the gap- 
less phase is stable: it forms stripes (as opposed to lines) 
in the phase diagram, bounded by ll s = ±6i' ± as 
shown in Fig. [TJ We have also computed the phase di- 
agram numerically from Eq. ([TO")) for the values t = 4, 
A = 2, /x = 0.5, which is also shown in Fig. [TJ the nu- 
merical data agree with our analytically obtained phase 
boundaries very well. 

V. CONCLUSIONS 

With the current interest in topological insulators and 
superconductors, the role of the block offdiagonal struc- 
ture [Eq. ([2])] in the literature so far has been in the topo- 
logical characterization of gapped phases of T-invariant 
superconductors. In this work we asked: what can be 
deduced from Eq. ([2]) regarding the nature of gapless 
phases of these systems? Focusing on systems in d — 2, 3, 
our central finding is that Eq. leads to the possibil- 
ity of d — 2 dimensional topologically stable Fermi sur- 
faces, characterized by an integer topological charge q. 
These Fermi surfaces can be envisioned as momentum 
space vortices in detD(k) as an "order parameter": in 
d = 2, they are points (in d = 3 they are lines) in momen- 
tum space where det-D(k) vanishes and their topological 
charge is calculated by the phase winding of det_D(k) 
around them. 

This result, combined with the momentum reversal re- 
lations in Eq. ([6]), leads to a "fcrmion quadrupling the- 
orem" : the number of Fermi points for d = 2 lattice 
models with a spinxEH structure is a multiple of four. 
Our analogous finding for d = 3 systems is that Fermi 
lines come in pairs. The most well known manifestation 
of the fermion quadrupling theorem is provided by the 
four Fermi points of a d-wave superconductor, and an 
example for a pair of Fermi lines is present in the non- 
centrosymmetric superconductor CePtsSiJi 

In Sec. [ITJ we have seen that our results provide a 
straightforward topological explanation for the robust- 
ness of the Fermi points of a d-wave superconductor 
against a weak spin-orbit coupling or against an admix- 
ture of a (time-reversal invariant) triplet pairing com- 
ponent. In addition, our considerations regarding two 
dimensional lattice systems can find applications similar 
to the work of Berg et aS^ who studied the perturba- 
tive stability of Fermi points in the presence of a (weak) 
coexisting order, with a periodicity commensurate with 
the lattice. For a T-invariant coexisting order, the au- 
thors of Ref. found that the Fermi points cannot be 
gapped out as long as the perturbation does not couple 
them to each other. With our topological considerations 
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at hand we can say a little bit more: the Fermi points 
are perturbatively unstable only if vortices are coupled 
to antivortices. In this case, the resulting Fermi point 
in the folded Brillouin zone has q = 0, which guaran- 
tees no topological stability; it can be gapped out by the 
perturbation. If Fermi points with like charges are cou- 
pled (for example a momentum reversed pair of Fermi 
points), the situation is different: in the folded Brillouin 
zone the resulting Fermi point has \q\ > 1. All what the 
small perturbation can do is to split it up to more Fermi 
points. Due to the quadrupling theorem, in addition, we 
know that the perturbation cannot remove only a single 
vortex-antivortex pair; it has to reduce the number of 
Fermi points by four. 

We have also found several implications of our findings 
to topological superconductors. For the surface of a topo- 
logical superconductor, we have related the net winding 
g no t to the bulk topological invariant v in the form of an 
index theorem, Eq. ([7]) . It accounts for the bulk topolog- 
ical properties in terms of a topological number for the 
boundary degrees of freedom, in perfect analogy to how 
the net number of chiral edge modes in a chiral supercon- 
ductor accounts for the Chern number— (the bulk topo- 
logical invariant for chiral superconductors). Eq. ([7]) tells 
us that v Fermi points on the surface are protected due 
to topological reasons, but it also shows that topology 
does not forbid to have more Fermi points in the surface 
theory. These extra Fermi points, if present, should be 
clearly visible, for example in the transport properties of 
the topological superconductor surface. 

Regarding the phase diagram for topological phase 
transitions, we have shown that generically, because of 



the presence of topologically stable Fermi surfaces, there 
is an extended gapless region in parameter space sepa- 
rating gapped phases with different v - gapless systems 
thus represent a separate, stable phase themselves. This 
should be contrasted to the case where there can be no 
topologically stable Fermi surfaces. In this case, gapless 
systems can only represent (unstable) phase boundaries. 

Finally, we note that our results related to topological 
superconductors can be extended to a bit broader context 
than the domains considered so far. First, as the symme- 
try classes in Eq. © contain all the possibilities 1 -^ in the 
Altland-Zirnbauer classificatio n 24 ! 25 for gapped three di- 
mensional systems with v € Z, our results can be seen as 
pertaining to " Z-topological materials" in d = 3, includ- 
ing chiral topological insulators 1 ^ in class AIII in addi- 
tion to superconductor systems. Second, as the key ingre- 
dient for the appearance of a stable gapless phase across 
a topological phase transition is the existence of topo- 
logically stable Fermi surfaces, such gapless phases are 
expected to appear also in two dimensional T-invariant 
topological superconductors. In these systems, the sta- 
bility of the gapless phase is provided by the topological 
stability of Fermi points. 
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